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On exact polytropic equilibria of self-gravitating gaseous 
and radiative systems: their application to molecular cloud 
condensation 
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ABSTRACT 

We propose a novel mathematical method to construct an exact polytropic sphere in 
self-gravitating hydrostatic equilibrium, improving the non-linear Poisson equation. 
The central boundary condition for the present equation requires a ratio of gas pressure 
to total one at the centre, which is uniquely identified by the whole mass and molecular 
weight of the system. The special solution derived from the Lane-Emden equation can 
be reproduced. This scheme is now available for modelling the molecular cloud cores in 
interstellar media. The mass-radius relation of the first core is found to be consistent 
with the recent results of radiation hydrodynamic simulations. 
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1 INTRODUCTION 

Fundamental research on the structure and evolution of 
stars, molecular clo uds, and galaxies is still very impor- 
tant iLongaiJ ^^). As a rule, a family of the Lane- 
Emden (LE) equations classified by the polytrope indices 
n is known to be the most powerful to ol to investisate the 
self-R ravitatiuK polytropic equilibria iBinnev &: Tremaind 
Il987l) . So far, many numerical and semi-analy tical 
works have been devoted to seeking solutions l| Liul 
ll996l:lRoxburgh fc Stockmanlll99gl:lGoenner fc Havasll200C 

;kl2 



trgri 

lHunteil200ll :fMedvede v fc Rvbickl200 J) , which include the 
non-singular function with n ^ oo for an isotherm al poly- 
tropic sphere le.g.. Nataraian fc Lvnden-Belll997^ . Remov- 
ing infinite divergence of the enclosed mass, some modified 
models are used to depict the globular clusters , intergalac- 
tic media, clusters of galaxies, and so forth (iKind Il962t 
fjones fc F orman 1984} . 

Prior to mod ern numerical methods 

JVandenBerg fc Bed Il985l: IVandenBerj |l98l, the pio- 
neering stellar model proposed by Eddington, which is 
based on the solution of the LE equation with the polytrope 
index of n = 3, could be useful to envisage stellar interiors, 
although the complexities such as transport processes and 
non-ideal equation of state (EOS) have been omitted. Re- 
garding the study of the interstellar medium (ISM), it seems 
that some classes of the LE equations shed light on dark 
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mole cular clouds and their condensations jMizuno et alJ 
I1994I) . However, in the general case that the pressure which 
supports the celestial objects is provided largely by the 
thermal motions of particles, we need to pay attention to 
the extensive use of the LE equation with a specified index. 
For heuristic ways, we explain this crucial point below. 

Let us suppose that the object sustains its own self- 
gravity due to the total pressure of P — pg -\- pi, where 
the partial pressures of gas and radiation field are de- 
scribed by the EOS of p^ = f3P = [kB/iiLmp)]{T/V) and 
Pr = (1 — /3)-P = l'^^^*! respectively. Here, mp is the proton 
mass, V = is the specific volume, a = 87r^feg/ {15c^h^) is 
the Stefan-Boltzmann constant, and ^ ^ 1 denotes the mean 
molecular weight. ^ For < /3 < 1, the relation between tem- 
perature and density is given by T = {[k-B/{iimp)]{3/a)[{l-~ 
f3) / /3]}^^^V~^^^ . Hence, we obtain the polytropic relation as 
follows; 



P(r) 



3 1-13 
a /34 



1/3 



V 



-4/3, 



r). 



(1) 



If we assume /3 = const dChandrasekhaJ ^G^), equa- 
tion (1) has the form of P = Kp^^^, where K = 
{[kB/{pmp)]'*'{3/a)[{l-/3)/l3*]y^'^ is the adiabatic constant. 
In this aspect, such a radiative polytrope leads to the LE 



1 For plasmas, 1/2 ^ /i ~ 2/(1 -)- 3X -|- 0.5Y) S 2, where the 
mass fractions contained in the unit mass of the medium are X g 
hydrogen, Y g helium, and Z = l — X — Yg heavy elements. 
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equation of index n = 3 for gravitational equilibria. Here- 
after, we refer to this scheme, namely, the Eddington model, 
as 'LE3', without further notice. Apparently, for gaseous ob- 
jects of /3 ~ 1, equation (1) does not match with the relation 
of P{r) oc p"'^^{r) involving n = 3/2 for a perfect monatomic 
gas. 

The possible scenario that recovers the physical con- 
sistency is to allow the spatial variation of polytrope in- 
dex. In equation (1) such changes can be achieved by in- 
voking the density dependence of /3 = /3[y(r)], and a triv- 
ial solution of dfi/fj, = (not shown). Of course, these 
are not mathematical artefacts. In the astronomical con- 
text, the present scheme is applicable to study, especially, 
the high-density adiabatic cores deeply immersed in the 
isothermal molecular cloud, which are so-called 'first cores' 
iMasunaea. Mivama fc Inutsukalll998^ . As yet, there has 
been no distinct detection of the first cores. We hope for 
future-planed observations of molecular cloud cores in the 
range of millimetre and submillimet re wavelengths with high 
angular resolution teasegawall200j) . 

In this paper, we present a self-consistent solution of a 
polytropic sphere for self-gravitating gaseous and radiative 
systems, going beyond the LE framework. In the following, 
we derive a generic equation for hydrostatic equilibria. In 
the present scheme, if we know the whole mass and the 
mean molecular weight of the system, then the pressure ra- 
tio of gas/total pressure at the centre, which contributes to 
the boundary condition, is identified. It is shown that the 
polytrope varies radially, modifying the density and pressure 
profiles derived from the LE3 model. In the limit that the 
gaseous pressure is asymptotically close to the total one, 
the solution exactly agrees with the LE solution of index 
n = 3/2. 



2 THE NONLINEAR POISSON EQUATION 
INCLUDING A SELF-CONSISTENT VARIABLE 
POLYTROPE 

We begin by defining the adiabatic exponent as Fi = 
d(lnP)/d(lnp). Introducing the total derivative of dP = 
d(Pg +Pr) = ~PgidV/V) + (pg +4pr)(dr/r), we obtain 



{Pg+4pr)^ - [pg - (Pg +Pr)ri]^ = 0. 



(2) 



For a quasi-static and adiabatic change, the thermodynam- 
ical principle requires 



dQ 
V 



-Pg + 12p, 



— + {pg + 4p,) — 



0, 



(3) 



where 7 = Cp/Cv is the specific-heat ratio. Combining equa- 
tion (2) with equation (3), we find the relation between Fi 
and (3 in the form of 



r^mV),'y]=f3{V) + 



(7-1) [4-3/3(^)1^ 
/3(1/) + 12 (7 - 1) [l-f3{V)Y 



(4) 



It should be noticed that equation (4) holds the asymptotic 
properties of Fi — > 4/3 and 7 for /3 — >■ and 1, respectively. 

On the other hand, taking the derivative of equation (1), 
we obtain dP/dV = -(P/3){(4- 3/3)/[/3(l - /3)]}(d/3/d1/) - 
(4/3) (P/y). This may be rewritten as 



dP 



4dV 



+ 



3/3 d/3 



3 V 3(1-/3) /3 



= 0. 



According to the relation of dP/P + TidV/V ■ 
we obtain the new relation of 



V 



d/3 _ 3/3(1 - /3) 



dV 



3/3 



ri(A7) 



(5) 

0, therefore, 
(6) 



which effectively changes the polytrope index of equation 
(1). 

Now we consider a hydrostatic equilibrium of a self- 
gravitating system, described by the spherically symmetric 
Poisson equation, namely, V"l> — r~^{d/dr){r^$) = AivGp, 
where $(r) — — p~^(dP/dr) and G is the gravitational con- 
stant. By using equations (1) and (6), the non-linear equa- 
tion can be cast to 



1 

4^ 





)'-] 


1/3 




1 a 






l-/3(p)" 







1 1/3 



Fi [/3(p),7]dp(r) 



dr 



p{r). (7) 



In order to normalize equation (7), we introduce p(^) = 
pis,)/ pc, where pc is the central mass density, and we prop- 
erly choose the dimensionless frame of ^ = r/a. Here, the 
characteristic scalelength a is given by 



a{pc,p) 



(47rG')3pi 



1/6 



(8) 



Then, we newly find a set of ordinary differential equations 
in the following form: 



1 dm(^) 



MO = ^ 



mo 



p2/3(e) dc 

1 dp(C) _ 



1 ~ m) 



no 
4 - 



-,1/3 

ri[/3(0,7], 

mo mo 



piO dC 3/3(C)[l-/3(C)]{Fi[/3(0,7] 



(9) 



Note that the enclosed mass defined by m(r) = 
Att p{r')r'^dr' is transformed into miO = 
—4'jTpca'^{pc, p)'rh{£,) ^ 0. In the special case of d/3/dp ~ 0, 
equation (9) reduces to the LE equation with an index of 
n ~ (7 — 1)^^ and 3 for /3 ~ 1 and 0, respectively. 



3 NUMERICAL SOLUTIONS OF THE 
GENERALIZED EQUATION 

To survey the inner structure of the polytropic sphere, we 
numerically integrate equation (9) for 7 = 5/3, and the 
boundary conditions of p{0\i^o = 1, '^(C)l5->o ~ 0, and 
/3(^)|5^o = /3c where /3c = {pg/P)c is a parameter given at 
the centre. It is noted that the solutions of p(^) and m(^) 
tend to monotonically decrease outwards, while /3(^) mono- 
tonically increases. At ^ = ^1, p(^) and drh{(,)/d(, vanish, 
indicating a well-defined radius and mass of the object, i.e. 



R = a(pc,p)Ci(/3c 
M = m[ei(/3c)] = 



-Anp^a (pc,^)m[^i(/3c)]. 



(10) 
(11) 
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Table 1. Values of gi, m, fi, and pc as a function of /3c 



_ 100 



(1 - /3)c 




m(^l) 




X 10 ^'^ g cm ^ 


0.01 


2.72 


-1.10 


1.11 


8.62 


0.05 


3.77 


-2.76 


1.75 


9.17 


0.1 


4.55 


-4.49 


2.24 


9.95 


0.3 


7.57 


-14.2 


3.98 


14.4 


0.5 


12.2 


-38.6 


6.57 


22.1 


0.7 


21.6 


-1.36 X 10^ 


12.3 


35.3 


0.9 


56.8 


-1.48 X 10=* 


40.7 


58.6 



For fj, = 2, see also Fig. U] (dotted curve). 

Table 1 lists the typical values of and 7n(^i) as a function 
of /3c. 



3.1 Inner structure of gaseous and radiative 
polytropic sphere 

Substituting equation (8) into equation (11), we ob- 
tain the explicit notation of the mean molecular weight: 
^ = [(4^G3Af2)-i(fcB/mp)4(3/a)]i/4{-m[Ci(/3c)]}'/^ Tak- 
ing the allowed parameter region into consideration, this can 
be expressed as 



/i(M,/?c) = 1.06 



Mq 

M 



1/2 



{-rfiMm'^' ^ 1, 



(12) 



which corresponds to the LE3 solution of 
/XLE3 = 4.25(M0/M)i/2(^ _ /3^)i/4/^^^ Further- 
more, using equations (8) and (10), we obtain 
Pc = {(47rGi?")-'[fcB/(M'np)]''(3/a)}i/2^?(^^)^ Making 

use of expression (12), this can be written as 
Pc(M,i?,/3c) = 4.73 X 10"^* 



— "if 



■ 100 au 



R 



{-m[Ci(/3c)l} 



g cm 



(13) 



Thus, for given M and /i, the fraction of gas pressure 
/3c is identified by equation (12). In addition, for given 
R, the central mass density is determined by equation 
(13), in contrast to the LE3 solution of pc le3 ~ 7.68 x 
10"^^(M/Mq)(100 au/Rf g cm"^, which is independent 
on /3c. These dependencies are summarized in Table 1. 

In Fig. we plot equations (12) and (13) as a func- 
tion of (pr/P)c = (1 - /3)c. It is found that, for M = Mq 
as an example, the inequality of equation (12) requires the 
fraction of radiation pressure at the centre to take its pa- 
rameter region of (1 — /3)c ^ 6.7 x 10"'^. ^ If we spec- 
ify the mean molecular weight, then the fraction can be 
fixed numerically. For fj, = 2 (hydrogen molecule), we ob- 
tain (1 — /3)c = 7.4 X 10~^. In such a gas dominant regime 
of (1 — /3)c <S 1, the dimensionless quantity ~ pc{R?' /M) 
does not largely depend on /3c, as seen in the figure. For 
R = 100 au, the numerical solution reveals the central den- 
sity of Pc = 8.6 x lO"'^^ g cm~^, which is about 90 per cent 

2 For plasmas, replace with 4.4 X IQ-* ^ {I - /3)c ^ 7.4 X 10"^, 
corresponding to the allowed parameter region ofl/2^/x^2 
(see footnote 1). 
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Figure 1. Central mass density of pc (bold solid curve) and 
Pc LE3 (solid line), both multiplied by (M0/M)(_R/1OO au)^, 
and mean molecular weight of fi (bold dashed curve) and /^le3 
(dashed curve), both multiplied by (M/Mq)^/^, as a function of 
the central pressure ratio of {pr/P)c- For an explanation, see the 
text. 



lower than that of pc le3 = 7.68 x 10~^^ g cm"''. For mas- 
sive objects of M Mq, the radiation dominant regime 
appears. In the limit of /3c — > 0, both pc and /x are asymp- 
totically close to Pc LE3 and ple3, respectively. This charac- 
teristic is a result of the behavior of Fi — > 4/3 for /3 — + in 
equation (4). 

In Fig. m for 0.01 ^ (1 - /3)c S 0.9, we show the ra- 
dial profiles of the normalized mass density p{r)/ pc, enclosed 
mass m{r)/M, total pressure P{r)/Pc, as well as the adia- 
batic exponent ri(r), comparing the LE solutions of index 
n = 3 and 3/2. Here the radius r = is normalized by 
R. It is shown in Fig. HI a) that the density profile gradually 
deviates from that of the LE3 solution as (1 — /3)c decreases. 
Concerning the case that the parameters of M and R are 
fixed, in the lower (1 — /3)c the central density decreases, be- 
cause the density varies sufficiently slowly from the centre 
to the envelope. In the limit of /3c — > 1, the present solu- 
tion shows an exact agreement with that of the LE equation 
with a polytrope index of n = 3/2. The radial profiles of the 
enclosed mass are shown in Fig. I^Jb). For higher (1 — /3)c 
including the LE3, the mass is more condensed in the cen- 
tral region. For example, the radius confining the half-mass 
indicates -Rm/2 le3 = 0.325-R. As (1 — /3)c decreases, this ra- 
dius gradually shifts outward, up to Ryi/2 lei. 5 = 0.521_R. 
It is confirmed that the gaseous polytropic effects are likely 
to flatten the density profile significantly. 

In Fig.l^Ic), we show the radial profiles of T\[l3{^),"f — 
5/3], defined by equation (4). As the argument /3(5) mono- 
tonically increases outwards, the adiabatic exponent also 
monotonically increases, yielding the self-consistent variable 
polytrope. This is one of the most important results in 
the present paper. For all /3c being considered, the expo- 
nent varies quite slowly near the central region, while for 
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Figure 2. Radial profiles of (a) normalized mass density p{r)/pc, 
(b) normalized enclosed mass m(r)/M, (c) adiabatic exponent 
ri(r), and (d) normalized total pressure P{r)/Pc, for (pr/P)c = 
0.01 (solid curves), 0.1 (long-dashed curves), 0.3 (short-dashed 
curves), 0.5 (dot-dashed curves), 0.7 (triple dot-dashed curves) 
and 0.9 (dotted curves). The figures have a common horizontal 
axis indicating the normalized radius of r/R. For comparison, 
the LE solutions with the polytrope index of n = 3 (LE3; bold 
solid curves/line) and n = 3/2 (LEI. 5; bold dotted curves/line) 
are also plotted. Note that in (a), (b) and (d) the solid curves 
for {pr/P)c = 0.01 overlap the bold dotted curves for the LEI. 5 
solution (see text). 



(1 — /3)c ~ 0.5, the variation around the envelope tends to 
be very steep. 

In Fig. I^Jd) , we show the radial profiles of total pres- 
sure of P(C) = [GMV(47ri?*)]{-C?(/3c)/m[6(/3c)]}'{[l - 
/3(C)]//3*(0}^^^p''''^(C), normalized by the central pressure 
Pc- It is noteworthy that, for the gas dominant regime of 
(1 — /3)c <^ 1, the pressure formula exhibits the density de- 
pendence approximated by P(^) oc p^^^{^) with Fi ~ 7, 
fully consistent with Fig. |5Jc). At ^ = 0, the pressure takes 
the peak value, to give 



Pc = 4.19 X 10" 



-Vf 



100 au 
R 



(l-/?c) 



f3t 



1/3 



ergs cm 



(14) 



Notice that, in contrast with the LE3 solution of Pc le3 = 
O.582(M/M0)^(1OO au/P)* ergs cm"^, equation (14) does 
depend upon /3c, but weakly for (1 — /3)c <^ 1. For M — Mq, 
/i = 2, and P = 100 au, we obtain the central pressure of 
Pc = 4.6 X 10~^ ergs cm"'^, which is an order of magnitude 
smaller than Pc le3. 

Invoking the EOS, the temperature profile can 
be described as T(C) = [3GAfV(47raP'')]i/''Ci(/3c) 
x{-m[Ci(/3c)]}-^/^{[l-/3(e)]//3(C)}'/'p'/'(e) (not shown 
in figure). At ^ = 0, it takes the peak value of 



6 (/3c 



(1 - /3c) 



/3c 



1/3 



K, 



(15) 



solution reads Tc le3 ~ 3.90 x 
au/P)(l - Pcf^'^ K; both having the 
Mq, h = 2, and P = ICQ au. 



whereas the LE3 
10^(M/Mq)1''2(100 
dependence of /3c. For M 

the central temperature of equation (15) is found to be Tc 
1.1 X 10^ K, which is again lower than Tc le3 = 2.0 X 10^ K. 
It is found that the decrease of the temperature is relatively 
small, within a factor of 2 for each /3c. 



3.2 Their application to the molecular cloud 

condensation: the critical radius of the 'first 



In the context of the study of molecular cloud condensa- 
tion in the ISM, the present quasi-stationary model is now 
available to provide an insight into the complicated dynam- 
icsjjPenston 1969; Larson, 196 9; Shu 1977; Saigo & Hariawal 
Il99a) . in particular, the formation of the first core (Ma- 
sunaga et al. 1998). For such an application, the effects of 
magnetic fields and turbulence might be taken into account, 
and the additional pressure 5P might be effectively included 
in equation (1) by repla cing n with n' = njS-o/i'i + /3p), 
where /3p = Pg/<5P > (iBludman fc Kennedy! [l99^ : e.g. 
for pure magnetic pressure, /3p = (87rpg)/P^. Moreover, for 
Tc ^ 2000 K, the dissociation of hydrogen molecules triggers 
gravitational contraction of the quasi-stationary first core. 
By invoking the scaling of Tc oc l//i' and equation (15), 
therefore, we find the relation between the mass and radius 
of the first core as follows: 



P ^ 13 



l + /3p 
/3p 



M 



6 (/3c) 



0.05Mq 

(1 - /3c) 



{-mKi(/3c)]}i/2 



/3c 



1/2 



1/3 



(16) 



for fi — 2. Note that the constraint of equation (12) gives 
the dependence of /3c = /3c(M, /i = 2) in equation (16). 
On the other hand, in the LE3 we obtain the mass-radius 
(M-P) relation of P ^ 4.6(M/O.O5M0)/3c au, where /3c = 
/3c(M,^LE3 =2). 

In Fig. 13 for /3p ^ 1, we plot the numerical solution 
of equation (16), that is, the allowable parameter region 
in the M—R plane. For convenience sake, we also plot the 
fractional radiation pressure at the centre as a function of 
the mass of the core: (pr/P)c = 1 — /3c(M, — 2). For 
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Mass of the first core M/M. 

Figure 3. Minimum possible radius of the first core of hydrogen 
molecular clouds (solid curve) and fractional radiation pressure 
at the centre (dotted curve) as a function of the mass of the core. 
The upper-left side of the solid curve (R > iJniin) corresponds 
to the stable region, while the opposite side {R < -Rmin) cor- 
responds to unstable region, concomitant with the dissociation 
of hydrogen molecules. For comparison, the M—R relation from 
radiation hydrodynamic simulations is indicated by the shaded 
inset (Masunaga et al. 1998). 

M S> Mq, corresponding to /3c ^ 1 for fixed /i (com- 
pare Table 1), the minimum core radius of equation (16) 
can be well approximated by that from the LE3: f?min ~ 
1.9xlO^[(l-h/3p)//3p](M/lOOM0)^/2 au. On the other hand, 
for M <^ Mq, we numerically find that the minimum pos- 
sible radius of equation (16) scales as 

and (1 - /3)c ^ 2.8 x 10-*[/3p/(l + /3p)]*(M/0.05Mq)^ The 
estimation of (17) reasonably supports the results of radia- 
tion hydrodynamic calculations by Masunaga et al. (1998). 
They show that the mass and radius of the first core are 
M ~ O.O5A/0 and ii ~ 5 au > iimin, respectively, and the 
results do not largely depend on the initial conditions of 
wide parameter ranges. If J? ^ Rmin, the core tends to be 
gravitationally unstable and evolves along the dynamically- 
contra cting track, to self-organize th e 'second core' as a pro- 
tostar iMasunaea fc Inutsukallioool) . 



4 CONCLUDING REMARKS 

In conclusion, we have developed a generic scheme to con- 
struct an exact polytropic sphere of self-gravitating gaseous 
and radiative medium. The particular results derived from 
the newly modified Poisson equation for hydrostatic equilib- 
ria are: 

(i) the numerical solutions show that for all cases, the 



central density and pressure are lower than those from the 
LE function with the index of n = 3; 

(ii) the adiabatic exponent monotonically increases radi- 
ally outwards; and 

(iii) in the special case that the central polytrope is 
gaseous closely, the present solution reproduces the prop- 
erties of the LE function with n — 3/2. 

Within this framework, the whole mass of the system is con- 
nected with the central density, temperature, and the mean 
molecular weight. 

For an application to modelling the molecular cloud 
condensation in the ISM, we have newly found the scal- 
ing law of the critical radius of the first core. The pre- 
liminary result is in consistent with that of the radia- 
tion hydrodynamic simulations. We expect that the ma- 
jor consequence can be also referred to, for example, mu- 
tatis mutand is, the study of the passive phase of protoplan- 
etary discs illon da fc Nakaeawal Il999l) . stellar modelling 
jBasu. Pinsonneauh fc BahcalJl200d) . and so on. 
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